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I  SO 


I.  INTRODUCTION 


In  order  to  investigate  the  scattering  effects  associated  with 
antennas  in  the  presence  of  ship  structures,  the  following  Fortran  IV 
computer  code  was  developed  for  (NOSC)  Naval  Ocean  Systems  Center, 

San  Diego,  California.  The  computer  code  is  used  to  compute  the  far  zone 
scattered  fields  for  antennas  radiating  in  the  near  zone  of  structures 
made  of  flat  plates.  In  its  present  form  this  code  simulates  a  ship 
structure  by  a  set  of  finite  flat  plates  forming  a  convex  structure  for 
which  the  scattering  from  one  flat  plate  to  another  is  negligible.  Using 
the  present  code,  one  can  treat  the  structure  by  a  single  flat  plate,  a 
rectangular  box,  a  rectangular  pyramid,  etc.  This  is  the  first  of  a  gener¬ 
ation  of  computer  codes,  and  it  is  limited  to  the  simpler  scattering 
mechanisms.  The  final  code  after  a  three  year  study  will  allow  one  to 
model  ship  structures  by  arbitrarily  placed  finite  flat  plates  and  finite 
elliptic  cylinders.  This  final  model  will  allow  for  a  simulation  of  ship 
structures  in  the  presence  of  UHF  and  above  antennas  with  all  significant 
mutual  interactions  taken  into  account. 

The  present  code  is  limited  to  one  structure  which  can  be  simulated 
by  as  many  as  14  plates.  This  is  based  on  the  array  dimensions  in  the 
code  and  not  a  limitation  of  the  theory.  Each  plate  can  consist  of  6 
corners;  however,  each  corner  must  lie  in  a  plane  or  the  computer 
code  will  abort.  The  definition  of  the  plates  is  made  by  first  setting 
up  a  fixed  cartesian  coordinate  system  relative  to  the  structure  under 
investigation.  The  plates  are,  then,  defined  by  the  location  of  the 
corners.  The  antenna  location  is,  also,  specified  in  the  same  co¬ 
ordinate  frame.  One  should  note  that  the  fixed  coordinate  system  should 
be  chosen  such  that  one  can  easily  define  the  structure.  The  program 
has  the  flexibility  to  handle  arbitrary  pattern  cuts  relative  to  this 
coordinate  system  as  is  discussed  later. 

The  antenna  presently  considered  in  the  computer  code  is  simu¬ 
lated  by  a  set  of  electric  or  magnetic  elemental  radiators.  There 
is  a  maximum  of  six  such  radiators  which  is  limited  by  the  computer 
code  dimension  and  not  the  theory.  Each  electric  or  magnetic  radiator 
has  a  cosine  distribution,  arbitrary  length,  arbitrary  magnitude  and 
phase,  and  arbitrary  orientation.  This  elemental  antenna  is  considered 
initially  but  can  be  easily  modified  in  that  the  code  is  modular  in  con¬ 
struction.  In  this  case,  the  SOURCE  subroutines  can  be  easily  exchanged 
with  another  antenna  pattern  subroutine. 

The  present  form  of  the  computer  code  is  not  large  in  terms  of  com¬ 
puter  storage  and  executes  a  pattern  in  short  order.  The  storage  is,  of 
course,  dependent  on  the  dimensions  which  might  vary;  however,  the  present 
code  requires  approximately  100  K  bytes.  It  will  run  a  pattern  cut 
of  360  points  for  a  flat  plate  structure  with  one  antenna  in  approximately 
10  seconds  on  a  CDC-6600  computer. 


The  limitations  associated  with  the  computer  code  results  from  the 
basic  nature  of  the  analysis.  The  solution  is  derived  using  the 
Geometrical  Theory  of  Diffraction  (GTD)  technique  which  is  a  high  fre¬ 
quency  approach.  In  terms  of  the  scattering  from  a  finite  flat  plate, 
this  means  that  each  plate  should  have  edges  at  least  a  wavelength  long. 
In  addition,  antenna  elements  should  not  get  closer  than  about  a  wave¬ 
length  to  any  edge.  In  some  cases,  the  previous  wavelength  limit  can 
be  reduced  to  a  quarter  wavelength. 

The  plot  package  necessary  to  generate  the  plots  that  are  shown 
to  illustrate  the  use  of  the  code  is  not  an  integral  part  of  the  com¬ 
puter  code.  This  results  from  the  variations  in  basic  plotting  routines 
from  one  computer  system  to  another.  These  variations  make  it  presently 
impractical  to  write  a  general  plot  package  to  be  used  with  the  computer 
code.  However,  a  polar  plot  routine  is  available  for  use  with  the 
code  which  in  most  cases  can  be  easily  modified  to  run  on  most  com¬ 
puter  systems. 


TABLE  I 


BLOCK  DIAGRAM  OF  THE  INPUT  DATA  ORGANIZATION 
FOR  THE  COMPUTER  CODE 


II.  DEFINITION  OF  INPUT  DATA 


As  stated  earlier,  the  input  data  is  defined  relative  to  a  fixed 
cartesian  coordinate  system  which  is  selected  by  the  user  to  make  the 
definition  of  the  structure  as  simple  as  possible.  In  this  system, 
all  linear  dimensions  are  given  in  terms  of  inches  and  all  angular 
dimensions  in  degrees.  The  angular  quantities  are  expressed  in  terms 
of  theta  and  phi,  the  usual  spherical  angles,  defined  relative  to 
the  operator's  fixed  coordinates. 

It  is  felt  that  the  maximum  usefulness  of  the  computer  code  can 
be  achieved  using  it  on  an  interactive  computer  system.  As  a  con¬ 
sequence,  all  input  data  are  defined  in  free  format  such  that  the 
operator  need  only  put  commas  between  the  various  inputs.  This  allows 
one  on  an  interactive  terminal  to  input  data  without  being  concerned 
with  the  fixed  field  length  associated  with  a  fixed  format. 

The  organization  of  the  input  data  is  illustrated  in  Table  1.  Note 
that  all  read  statements  are  made  on  unit  number  5,  i.e.,  READ  (5,-), 
where  the  symbol  refers  to  free  format.  In  all  the  following 
discussions  associated  with  logical  variables  a  "T"  will  imply  true, 
and  an  "F"  will  imply  false.  The  complete  words  true  arid  false  need 
not  be  input  in  that  most  compilers  just  consider  the  first  character  in 
determining  the  state  of  the  logical  variable.  The  following  list 
defines  in  detail  the  function  associated  with  each  of  the  input 
variables : 


1. 


READ:  LDEBUG,  LTEST ,  LOUT 


a)  LDEBUG:  This  is  a  logical  variable  defined  by  T  or  F. 

It  is  used  to  debug  the  progran  if  errors  are 
suspected  within  the  program.  If  set  true,  the 
program  prints  out  data  on  unit  n 6  associated 
with  each  of  its  internal  operations.  These 
data  can,  then,  be  compared  with  previous  data 
which  are  known  to  be  correct.  It  is,  also,  used 
to  insure  initial  operation  of  the  code.  Only 
one  pattern  angle  is  considered. 

(normally  set  false) 

b)  LTEST:  This  is  a  logical  variable  defined  by  T  or  F. 

It  is  used  to  test  the  input/output  associated 
with  each  subroutine.  The  data  written  out  on 
unit  #6  are  associated  with  the  data  in  the 
window  of  the  subroutine.  They  are  written 
out  each  time  the  subroutine  is  called.  It  is, 
also,  used  to  insure  initial  operation  of  the  code. 
Only  one  pattern  angle  is  considered. 

(normally  set  false) 

c)  LOUT:  This  is  a  logical  variable  defined  by  T  or  F. 

It  is  used  to  output  data  on  unit  #6  associated 
with  the  main  program.  It  too  is  used  to  initially 
insure  proper  operation.  It  can  be  used  to  examine 
the  various  components  of  the  pattern. 

(normally  set  false) 

2.  READ:  LSLOPE,  LCORNR 

a)  LSLOPE:  This  is  a  logical  variable  defined  by  T  or  F. 

It  is  used  to  tell  the  code  whether  or  not  slope 
diffraction  is  desired  during  the  computation, 
(normally  set  true) 

b)  LCORNR:  This  is  a  logical  variable  defined  by  T  or  F. 

It  is  used  to  tell  the  code  whether  or  not 
corner  diffraction  is  desired  during  the  com¬ 
putation. 

(normally  set  true) 

3.  READ:  THC ,  PHC,  THP 

This  set  of  data  is  associated  with  the  conical  pattern  desired 
during  execution  of  the  program.  The  pattern  axis  is  defined  by  the 
spherical  angles  (THC,  PHC)  as  illustrated  in  Figure  1.  These  angles 
define  a  radial  vector  direction  which  points  in  the  direction  of  the 


pattern  axis  of  rotation.  These  angles  actually  set-up  a  new  coordi¬ 
nate  system  in  relation  to  the  original  fixed  coordinates.  The  new 
cartesian  coordinates  defined  by  the  subscript  "p"  are  found  by  first 
rotating  about  the  z-axis  the  anqle  PMC  and,  then,  about  the  y-axis  the 
angle  THC.  The  pattern  is,  then,  taken  in  the  "p"  coordinate  system 
in  terms  of  spherical  angles.  The  theta  angle  of  the  pattern  taken 
about  the  Zn-axis  is  defined  by  TMP.  The  phi  anqle  is  defined  by  the 
next  read  statement.  In  the  present  form  the  program  will,  then, 
compute  any  conical  pattern  in  that  TMP  is  used  as  the  conical  pattern 
angle  about  the  Zp-axis  for  the  complete  pattern  calculation. 

a)  THC,  PHC:  These  are  real  variables.  They  are  input  in 

degrees  and  define  the  axis  of  rotation  about 
which  a  conical  pattern  will  be  computed. 

b)  THP:  This  is  a  real  variable.  It  is  input  in  degrees 

and  used  to  define  the  conical  anqle  about  the 
axis  of  rotation  for  the  desired  pattern. 

4.  READ:  IB,  IE,  IS 

a)  IB,  IE,  IS:  These  are  integer  variables  used  to  define 

angles  in  degrees.  They  are,  respectively,  the 
beginning,  ending,  and  incremental  values  of  the 
phi  pattern  angle. 

As  a  result  of  the  input  given  by  the  two  previous  read  statements, 
the  operator  has  completely  defined  the  desired  conical  pattern  to  be 
computed  during  execution  of  the  proqram. 

5.  READ:  FRQG 

a)  FRQG:  This  is  a  real  variable  which  is  used  to  define  the 
frequency  in  gigahertz. 


6.  READ:  MPX 

a)  MPX:  This  is  an  integer  variable  which  defines  the  maximum 
number  of  plates  to  be  used  in  the  present  simulation 
of  the  structure.  Presently,  1  1  MPX  1  14. 

7.  READ:  (flEP(MP)  ,  MP  =  1  ,  MPX) 

a)  MEP  (MP):  This  is  a  dimensioned  integer  variable.  It  is 
used  to  define  the  number  of  corners  (or  edges) 
on  the  MPth  plate.  They  are  defined  on  the  same 
data  line  with  commas  between  integer  values. 
Presently,  1  _  MEP  (MP)  <_  6  with  1  MP  •_  14. 
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8.  READ:  (X(MP,ME,N),  N=l,3) 

As  stated  earlier  the  locations  of  the  corners  of  the  flat  plates 
are  input  in  tercis  of  the  x,  y,  z  coordinates  in  the  fixed  cartesian 
coordinate  system. 

a)  X(MP,ME,N):  This  is  a  triply  dimensioned  real  variable.  It 

is  used  to  specify  the  location  of  the  ME^h  corner 
of  the  MP^  plate  in  inches.  It  is  input  on  a 
single  line  with  the  real  numbers  being  the  x,y,z 
coordinates  of  the  corner  which  correspond  to 
N=l,2,3,  respectively,  in  the  array.  For  example, 
the  array  will  contain  the  following  for  plate  #1 
and  corner  #2  located  at  x=2.,  y=4.,  z=6.: 

X ( 1 ,2,1 )  =  2. 

X( 1 ,2 ,2 )  =  4. 

X( 1,2,3)  =  6. 

This  data  is  input  as:  2.,  4.,  6. 

Considering  the  flat  plate  structure  given  in  Figure  1,  the  input 
data  is  given  by 


-I 

"■I 


1.,  1.,  0. 
-1 . ,  1 . ,  0. 
-1.,  -1.,  0. 
1 . ,  -1 . ,  o. 


:  corner  #1 
:  corner  #2 
:  corner  #3 
:  corner  #4 


plate  #1 


Presently :  1  <_  MP  <  14 

1  <_  ME  1  6 

1  <  N  <  3 


(See  Section  III  for  further  details  on  defining  the  corner  points.) 
9.  READ:  MSX,  LSOR 

a)  MSX:  This  is  an  integer  variable  which  defines  the  maximum 

number  of  elemental  radiators  to  be  considered  during 
execution  of  the  program.  Presently,  1  ^  MSX  6. 

b)  LSOR:  This  is  a  logical  variable  which  is  defined  by 

T  or  F.  It  is  used  to  specify  whether  or  not  the 
operator  wants  simply  the  antenna  pattern  alone 
without  the  scattering  by  the  plate  structure. 


10.  REAO:  (XSS  (MS,N) ,  N=1 ,3) 


a)  XSS  (MS,N):  This  is  a  doubly  dimensioned  real  array  which 
is  used  to  define  the  x,y,z  location  of  the 
MSth  elemental  radiator  in  the  fixed  cartesian 
coordinate  system.  Again,  a  single  line  of  data 
contains  the  x,y,z  (N=l  ,2,3)  locations. 

Presently,  1  <  MS  <  6 
1  1  N  1  3 

(See  Section  III  for  further  details  on  defining  the  source  point.) 

11.  READ:  I  MS  (MS)  ,  HS  (MS),  THO  (MS),  PIIO  (MS) 

a)  IMS  (MS);  This  is  an  integer  array  which  is  used  to  define 

whether  the  MS^h  source  is  an  electric  or  magnetic 
elemental  radiator. 

IMS  (MS)  =04  electric 
IMS  (MS)  =  1  4  magnetic. 

b)  HS  (MS):  This  is  a  real  array  which  is  used  to  input  the 

length  of  the  MS*^  element  in  electrical  wavelengths. 

It  is  not  input  in  inches. 

c)  THO  (MS),  PHO  (MS):  These  are  real  arrays  which  are  used 

to  define  the  orientation  of  the  MSth 
element  in  the  fixed  cartesian  coordinate 
system.  The  THO,  PHO  angles  are  in  degrees 
and  define  a  radial  direction  which  is 
parallel  to  the  MS^  element  current 
flow.  Again,  THO  is  the  theta  angle  and 
PHO  is  the  phi  angle  in  a  normal  spherical 
coordinate  system. 

Presently,  1  <_  MS  6. 

12:  READ:  WM  (MS),  WP  (MS) 

a)  WM  (MS),  WP  (MS):  These  are  real  dimensioned  arrays  used 

to  define  the  excitation  associated  with 
the  MS^  element.  The  magnitude  is  given 
by  WM  and  the  phase  in  degrees  by  WP. 
Presently,  1  ^  MS  ^  6. 

This  concludes  the  definition  of  all  the  input  parameters  to  the 
program.  The  program  would,  then,  run  the  desired  data  and  output  the 
results  on  unit  #6.  However  as  with  any  sophisticated  program,  the 
definition  of  the  input  data  is  not  sufficient  for  one  to  fully  understand 
the  operation  of  the  code.  In  order  to  overcome  this  difficulty  the  next 
section  discusses  how  the  input  data  are  interpreted  and  used  in  the 
program. 
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III.  INTERPRETATION  OF  INPUT  DATA 


This  computer  code  is  written  to  require  a  minimum  amount  of  user 
information  such  that  the  burden  associated  with  a  complex  geometry 
will  be  organized  internal  to  the  computer  code.  For  example,  the 
operator  need  not  instruct  the  code  that  two  plates  are  attached  to 
form  a  convex  structure.  The  code  flags  this  situation  by  recognizing 
that  two  plates  have  a  common  set  of  corners  (i.e.,  a  common  edge). 

So  if  the  operator  wishes  to  attach  two  plates  together  he  needs  only 
define  the  two  plates  as  though  they  were  isolated.  However,  the 
two  plates  will  have  two  identical  corners.  All  the  geometry  information 
associated  with  plates  with  common  edges  is  then  generated  by  the  code. 

As  discussed  below  the  computer  code  has  two  subtle  interpretations 
which  must  be  kept  in  mind  while  definino  the  geometry.  The  first 
subtlety  is  associated  with  sirulatinq  convex  structures  in  which  two 
or  more  plates  are  used  in  the  computation.  This  problem  results  in 
that  each  plate  has  two  sides.  If  the  plates  are  used  to  simulate  a 
closed  or  semi-closed  structure,  then  possibly  only  one  side  of  the  plate 
will  be  illuminated  by  the  antenna.  Consequently,  the  operator  must  define 
the  data  in  such  a  way  that  the  code  can  infer  which  side  of  the  plate  is 
illuminated  by  the  antenna.  This  is  accomplished  by  defining  the  plate 
according  to  the  right-hand  rule.  As  one's  fingers  of  the  right  hand 
follow  the  edges  of  the  plate  around  in  the  order  of  their  definition,  his 
thumb  should  point  toward  the  illuminated  side  of  the  plate.  To  illustrate 
this  constraint  associated  with  data  format,  let  us  consider  the  definition 
of  a  rectangular  box.  In  this  case,  all  the  plates  of  the  box  must  be 
specified  such  that  they  satisfy  the  right-hand  rule  with  the  thumb 
pointing  outward  as  illustrated  in  Figure  2.  If  this  rule  were  not 
satisfied  for  a  given  plate,  then  the  code  would  assume  that  the 
antenna  is  within  the  box  as  far  as  the  scattering  from  that  plate 
is  concerned. 

The  second  situation  which  must  be  kept  in  mind  is  associated 
with  antenna  elements  mounted  on  a  plate.  First,  a  plate-mounted 
antenna  must  be  positioned  close  (approximately  10"°  wavelengths  or  less) 
but  not  exactly  on  the  plate.  If  the  antenna  is  positioned  exactly  on  the 
plate,  the  code  can  not  decide  which  side  of  the  plate  the  antenna  is 
mounted  on  for  the  calculation.  For  example,  a  monopole  is  mounted  on 
one  side  of  a  plate;  yet,  the  code  must  recognize  which  side.  This  point 
will  be  further  discussed  in  the  following  section.  There  is  another 
point  associated  with  antennas  mounted  on  perfectly  conducting  structures. 
If  a  plate-mounted  monopole  is  considered  'n  the  computation,  then  one 
should  input  the  equivalent  dipole  length  and  not  the  monopole  length. 

The  program  automatically  handles  the  half  dipole  modes  associated  with 
the  monopole.  Finally,  if  a  pi ate-mounted  slot  antenna  is  considered  in 
the  computation,  the  code  automatically  doubles  the  radiation  due  to 
the  image  effect  of  the  perfectly  conducting  ground  plane. 

The  following  section  uses  a  set  of  sample  problems  to  illustrate 
how  the  operator  can  realize  the  versatility  of  the  code  and  still 
satisfy  the  few  constraints  associated  with  the  input  data  format. 


Figure  2.  Data  format  used  to  define 


a  box  structure. 


i 


IV.  APPLICATION  or  CODE  TO  SFVLl'.AL  EXAMPLES 


The  followinq  exarip  es  arp  uspd  to  illustrate  the  various  features 
of  the  conputer  code.  The  results  are  presented  in  terns  of  polar  plots 
which  can  be  generated  usinq  the  polar  plot  routine  qiven  in  Appendix  I. 
Note  that  the  patterns  are  plotted  in  decibels  with  each  division  beinq 
10  dB.  The  routine  in  Appendix  I  does  not  nake  any  reference  to  a 
symbol  subroutine,  in  that  they  are  very  much  machine  dependent.  Thus, 
labeling  is  not  included  on  the  plots. 

1.  Consider  the  pattern  of  an  electric  dipole  in  the  presence  of  a 
finite  ground  plane  as  shown  in  Fiqure  3.  The  input  data  is  given  by: 

F.F.F 

T.T 

0. ,0. ,90. 

0,360,1 

10.43 

1 

4 

0. ,6.  ,6.  :  corner  ?1  \ 

0.  ,-6.  ,6. 

0. ,-6 . ,-6. 

0. ,6. ,-6. 

l.F 

6.75,0. ,0. 

0..1 ,0. ,0. 

1.  ,0. 

The  Egp  pattern  is  plotted  in  Finure  4.  The  F$p  pattern  is  not 
plotted  in  that  it  is  of  negligible  value. 


:  corner 
:  corner  *3 
:  corner  #4 


/  Plate  *1 

i 


4.  Eep  radiation  pattern  for  a  short  electric  dipole  located 
above  a  square  plate  at  a  frequency  of  10.43  GHz. 


2.  Consider  the  pattern  of  a  magnetic  dipole  in  the  presence  of 
finite  ground  plane  as  shown  in  Figure  3.  The  input  data  is  given  by 

F.F.F 

T,T 

0. ,0. ,90. 

0,360,1 

10.43 

1 

4 

0. ,6.  ,6. 

0. ,-6. ,6. 

0. ,-6. ,-6. 

0. ,6. ,-6. 

1  »F 

6.75,0. ,0. 

1 ,0.1,0. ,0. 

1 .  ,0. 


The  E<pp  pattern  is  plotted  in  Figure  5.  The  Eep  pattern  is  not 
plotted  in  that  it  is  of  negligible  value. 


h^p 


Figure  5.  E^p  radiation  pattern  for  a  short  magnetic  dipole  located 
above  a  square  plate  at  a  frequency  of  10.43  GHz. 
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3.  Consider  the  H-plane  pattern  of  an  electric  dipole  in  the 
presence  of  a  roof-top  type  structure  as  shown  in  Figure  6.  The 
input  data  is  given  by: 


F.F.F 

T.T 

90.  ,0. ,90. 
0.360,1 


4,4 

6.496,0. ,0. 
6.496,3.34,-3.34 
-6.496,3.34,-3.34 
-6.496,0. ,0. 
6.496,0. ,0. 
-6.496,0. ,0. 
-6.496,-3.34,-3.34 
6.496,-3.34,-3.34 
1  »F 

0. ,0. ,3.3 
0, .5 ,90. ,0. 
l.,0. 


corner  #1 
corner  #2 
corner  #3 
corner  #4 
corner  H] 
corner  HZ 
corner  #3 
corner  #4 


Plate  #1 


Plate  #2 


The  E0p  polar  pattern  is  plotted  in  Figure  7.  The  same  patterr 
is  plotted  in  rectangular  form  in  Figure  8  where  it  is  compared  wit! 
a  measured  pattern. 


Figure  7.  E©p  radiation  pattern  for  a  A/2  electric  dipole  located 
above  a  roof-top  type  structure  at  a  frequency  of  4  GHz. 


4.  Consider  an  off-principal  plane  pattern  of  an  electric  dipole 
in  the  presence  of  a  roof-top  type  structure  as  shown  in  Fiqure  9.  Th 
input  data  is  given  by: 

f,f,f 

T,T 

30. ,0. ,90. 

0,360,1 

4. 


4,4 

6.496,0. ,0. 
6.496,3.34,-3.34 
-6.496,3.34,-3.34 
-6.496,0. ,0. 
b. 496,0. ,0. 
-6.496,0. ,0. 
-6.496,-3.34,-3.34 
6.496,-3.34,-3.34 
1,F 

0.  ,0. ,3.3 
0,.5,90. ,0. 
l.,0. 


The  Eep  and  E^,n  polar  patterns  are  plotted  in  Figure  10.  The 
sane  patterns  are  plotted  in  rectanqular  form  in  Figure  11  where  they 
are  compared  with  measured  results. 


Figure  10a.  Eep  radiation  pattern  for  a  \/2  electric  dipole  located 
above  a  roof-top  type  structure  at  a  frequency  of  4  GHz. 


5.  Consider  a  principal  plane  pattern  of  an  electric  dipole  in 
the  presence  of  a  pyramid  structure  as  shown  in  Fiqure  12.  The  input 
data  is  given  by: 

f.f.f 

T.T 

90..  90.. 90. 

0,360,1 

11.81 

4 

3, 3, 3, 3 

6. . 0. .0. 


0. 

»6. , 

-6. 

-6 

.,0. 

,0. 

-6 

.  ,0. 

,0. 

0. 

,6., 

-6. 

0. 

,-6. 

,-6. 

-6 

.  ,0. 

,0. 

0. 

>-6. 

,-6. 

6. 

,0., 

0. 

6. 

,  0 . , 

0. 

0. 

-6., 

-6. 

0. 

,  6 . , 

-6. 

1. 

F 

0. 

,0., 

3. 

0, 

0.5, 

90., 

1. 

,0. 

The  E^p  pattern  is  plotted  in  Figure  13.  The  Eqp  pattern  is  not 
plotted  in  that  it  is  of  negligible  value. 


PLATE  #4 


Figure  12.  A  A/2  dipole  in  the  presence 
of  a  pyramid  structure. 


b.  Consider  the  pattern  of  a  a/4  nonopole  Mounted  in  the  center 
of  a  square  ground  plane  as  shown  in  Fiaure  14.  The  input  data  is 
given  by: 

f.f.f 

T»T 

90..  90.. 90. 

0,360,1 

10. 

1 

4 

6. . 6. .0. 

-6. ,6. ,0. 

-6. ,-6. ,0. 

6. ,-6. ,0. 
l.F 

0.,0.,1.E-10 

0, . 5 ,0. ,0. 

l.,0. 

The  E^p  pattern  is  plotted  in  Figure  15.  The  E@p  pattern  is  not 
plotted  in  that  it  is  of  negligible  value. 
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Figure  14.  A  A/4  monopole  mounted  on  a 
square  ground  plane. 


7.  Consider  the  pattern  of  a  A/2  slot  antenna  mounted  in  the  center 
of  a  square  ground  plane  as  shown  in  Figure  16.  The  input  data  is 
given  by: 


F.F.F 

T.T 

90.  ,135.  ,90. 
0,360,1 
10. 

1 

4 

6.  ,6. ,0. 

“6. ,6. ,0. 

-6. ,-6. ,0. 

6* ,-6. ,0. 

1  »F 

0. ,0. ,1 . E- 1 0 
1  ,.5,90. ,90. 
l.,0. 


The  Eep  and  E^p  patterns  are  plotted  in  Figure  17.  Note  this 
pattern  is  taken  directly  across  the  corner  of  the  plate;  yet,  both 
patterns  are  continuous. 
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Figure  17a.  E@p  radiation  pattern  of  A/2  slot  antenna  mounted  on  a 
square  ground  plane  at  a  frequency  of  10  GHz. 
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8.  Consider  the  pattern  of  a  A/2  electric  dipole  located  near 
a  rectangular  box  as  shown  in  Figure  18.  The  input  data  is  given  by: 

F.F.F 

T.T 

90. . 135. .90. 

0,360,1 

4. 

6 

4, 4,4, 4, 4, 4 
0.  ,0. ,0. 

6. . 0. .0. 

6. ,6. ,0. 

0. ,6. ,0. 

0. ,6. ,0. 

6. ,6. ,0. 

6, ,6* ,-6. 

0. ,6. ,-6. 

0. ,0. ,-6. 

0. ,6. ,-6. 

6. . 6.  i*6 . 

6. . 0. ,-6. 

6. ,0. ,0. 

0. ,0. ,0. 

0. ,0.  ,-6. 

6. ,0.  ,-6 . 

6..6..0. 

6. ,0. ,0. 

6. ,0. ,-6. 

6. ,6. ,-6. 

0.  ,0. ,0. 

0. ,6. ,0. 

0. ,6. ,-6. 

0. ,0. ,-6 . 
l.F 

-1.575,-1.575,2.23 
0, .5,90. ,-45. 

l.,0. 


The  Eep  polar  radiation  pattern  is  plotted  in  Figure  19.  Note 
that  the  E^p  is  of  negligible  value.  This  same  pattern  is  plotted 
in  rectangular  form  in  Figure  20,  where  it  is  compared  with  a  measured 
result. 
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>  36  72  108  144  180  216  252  288  324 

0  (DEGREES) 


Figure  20.  Comparison  of  measured  and  calculated  E@p  radiation 
pattern  for  a  dipole  near  a  box  in  the  indicated 


9.  Consider  the  various  radiation  patterns  of  an  array  of  A/2 
electric  dipoles  situated  1/4"  above  a  square  ground  plane  as  shown 
in  Figure  21 . 

a)  Input  data  for  6  element  array  scanned  at  broadside: 

F.F.F 

T.T 

-90. ,0. ,90. 

0,360,1 

10. 

1 

4 

6. ,6. ,0. 

-6. ,6. ,0. 

—6. ,— 6. ,0. 

6. ,— 6. ,0. 

6,F 

0. ,-l .25,0.25 
0,0.5,90. ,0. 
l.,0. 

0. ,-.75, .25 
0, .5,90. ,0. 

l.,0. 

0. ,-.25, .25 
0,.5,90. ,0. 

l.,0. 

0. ,.25, .25 
0,.5,90. ,0. 

l.,0. 

0. ,.75, .25 
0, .5 ,90. ,0. 

l.,0. 

0., 1.25, .25 
0,0.5,90. ,0. 

l.,0. 

The  Eep  radiation  pattern  is  plotted  in  Figure  22.  Note  that  the 
ground  plane  has  little  effect  for  the  broadside  case. 


Figure  21.  Six  element  A/2  dipole  array  with  uniform  amplitude 
and  phase  scanned  above  a  square  ground  plane. 


\-^4>P 


b)  Input  data  for  G  element  array  scanned  ^36  from  broadside 


F.F.F 

T.T 

-90. ,0. ,90. 
0,360,1 
10. 

1 

4 

6. ,6. ,0. 

— 6.  ,6. ,0. 

-6. ,-6. ,0. 

6. ,— 6. ,0. 

6,F 

0. ,-1.25,0.25 
0,0.5,90. ,0. 

1  .  ,225. 

0.  ,-.75, .25 
0,.5,90. ,0. 

1 .  ,135. 

0.  ,-.25, .25 
0,.  5,90. ,0. 

1 .  ,45. 

0. ,  .25,  .25 
0,. 5,90. ,0. 

1 .  ,-45. 

0.  ,.75, .25 
0,.5,90. ,0. 

1 .  ,-135. 

0., 1.25, .25 
0,0.5,90. ,0. 

1 . ,-225. 


The  Egp  radiation  pattern  is  plotted  in  Figure  23. 


Eq  p  radia 
electric  < 
plane  sea 
10  GHz. 


pattern 
les  situa 
^36°  o 


c)  Input  data  for  6  element  array  scanned  V2&  from  broadside 


F.F.F 

T.T 

-90.  ,0.  ,90. 
0,360,1 
10. 

1 

4 

6. ,6. ,0. 

— 6. ,6. ,0. 

—6 . , —6 . ,0 . 

6. ,-6. ,0. 

6,F 

0.  ,-1.25,0.25 
0,0.5,90. ,0. 


1. 

,389.71 

0. 

,-.75, 

.25 

0, 

.5,90. 

,0. 

1. 

,233.83 

0. 

,-•25, 

.25 

0, 

.5,90. 

,0. 

1. 

,77.94 

0. 

,.25,. 

25 

0, 

.5,90. 

,0. 

1. 

,-77.94 

0. 

75 

25 

0, 

.5,90. 

,0. 

1. 

,-233. 

83 

0. 

,1.25, 

.25 

0, 

0.5,90 

.,0 

1. 

,-389. 

71 

The  Eep  radiation  pattern  is  plotted  in  Figure  24. 


V. 


PROGRAM  OUTPUT 


The  basic  output  from  the  computer  code  is  a  line  printer  listing 
of  the  results.  Recall  that  the  results  of  the  program  are  the  Eep 
and  E^p  radiation  pattern  values.  In  order  to  again  describe  these 
pattern  components,  let  us  consider  the  first  few  examples  of  the 
previous  section.  The  computer  code  allows  for  a  rotation  of  coordinates 
in  terms  of  the  pattern  calculation.  This  is  input  in  terms  of  the 
spherical  angles  THC,  PHC.  Using  the  geometry  of  examples  *1  and  #2  as 
given  in  Figure  3,  THC  =  PHC  =  0°,  such  that  the  pattern  coordinates  are 
aligned  with  the  structural  coordinates.  Thus.  Eep  wou)d  be  the^normgil 
theta  component  of  the  field  (i.e.^Eep  =  £  ;  ep  =  f  •  e,  since  6p  =  e 
in  this  case).  Likewise,  E^p  =  £  4p  =  f  •  $  in  this  case.  Note  that 
£  is  the  total  radiated  electric  field. 

Now  let  us  suppose  that  one  defines  the  same  structure  as  above 
using  the  coordinate  system  shown  in  Figure  25.  In  order  to  take  the 
same  pattern  cut  as  in  the  previous  case,  one  must  set  THC  =  90°  and 
PHC  =  180°.  Then  with  THP  =  90°,  one  will  obtain  the  same  pattern  as 
given  previously;  however  in  this  case,  the  pattern  data  will  be 
inverted  in  terms  of  the  plot  as  shown  in  Figure  26.  Again,  the  radiation 
pattern  is  output  in  terms  of  Eep  and  E^p  as  given  by 


where  0p  and  <j>p  are  the  normal  spherical  angles  in  the  xp,  yp,  zp  co¬ 
ordinate  system. 

In  order  to  further  understand  this  concept  associated  with  the 
rotated  pattern  coordinates,  the  pattern  coordinate  system  (xp,  yp,  zp) 
is  illustrated  in  each  of  the  figures  given  in  the  previous  example 
section. 

In  that  one  is  not  able  to  deliver  standard  plot  routines  from 
one  computer  system  to  another,  our  plot  package  is  not  included 
as  an  integral  part  of  this  computer  code.  However,  a  simple  polar 
plot  routine  is  given  in  Appendix  I  which  can  be  used  to  duplicate 
the  polar  graphs  given  in  the  previous  section.  This  plot  routine  is 
included  as  part  of  the  code  such  that  one  must  be  mindful  of  its 
inclusion.  It  is  possible  that  the  program  will  not  compile  if  a 
"PLOT"  subroutine  is  not  available. 
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Figure  25.  Example  of  the  pattern  coordinates  using  a  different 
fixed  coordinate  system  for  the  structure  of 
Figure  3. 


I 
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o  c.  cj  a  o  o  i  i  o  c  c ;  noo  o  o  o  r. 


APPENDIX  I 


The  following  is  a  listing  of  a  polar  plot  subroutine  which  can  be 
used  to  generate  the  polar  patterns  given  in  Section  IV.  Note  that 
this  code  refers  to  two  subroutines  "PLOTS"  and  "PLOT"  which  must  be 
added  by  the  system  or  operator.  The  definitions  of  these  routines 
are  given  in  the  comments  associated  with  the  code. 


SUBROUTINE  P()LPLT<  ET,  RP ,  I  PLT,  I  PHS ,  I  DM  ) 

i  i 

!  j  THIS  ROUTINE  IS  USED  TO  PLOT  THE  RESULTS  II’  TERNS  OF  A 
!!  POI  AR  PLOT.  TUF  C/II  TO  SUBROUTINE  "PLOTS"  IS  USED  TO  INITIALIZE 
!!  THE  PLOTTER.  VJF  CALLS  TO  THE  SUBROUTINE  "PLOT"  ARE  USED  TO 
! !  DRAW  THE  DATA  AS  FOLLOWS* 

;  i 

! !  CALL  PLOT ( X , V , N) 

i  ; 

!!  X ,  Y=COORDI  NATES  OF  THE  NFL  PLOT  P()!>-T. 

i  i 

ii  N=?  pcy  15  on WN  MOVING  TO  THE  MEW  POP'T. 

!!  N=3  PEN  IS  U°  MOVING  TO  TUE  MFW  POPT. 

!!  M=9H9  E"FFER  USED  TO  STORE  PLOT  r'AT*  IS  EMPTIED  TO  PLOTTER 

i  ; 

!  !  N<0  I'/OLIFS  ORIGIN  SHIFT  TO  THE  MEW  POINT. 

!  !  N>0  IMPLIES  NO  ORIGIN  SUIFT  AFT r*R  MOVING  TO  NEW  POINT. 

;  i 

COMPLEX  ETC  ID?') 

DIMENSION  P’UH(ICC) 

DATA  PI  ,  TPI  .DPR/S.  I  A  I  5926b.  6.  2><3  I  S55,  57  .295  7795H/ 

E  MX  =  0 . 

DO  101  I P=0,360, I PHS 
I-IP+1 

EM=CABS(HT( I  ) ) 

IFCEM.GT.EMX)  EMX=EM 
101  CONTINUE' 

CALL  PLOTSC II UE, 1 00.3 ) 

CALL  PLOT (4. 25,5.i>, -3) 

C  ***  POLAR  GRID  *** 

DO  110  1  =  1,/- 
R0=RP*I/4 . 

CALL  PLOTCRG.O.  ,3  ) 

DO  1 10  J=0, 36  0,2 
A  NG=J /DPR 
>  X=RC*COS  (  A  \'(  ) 

YY~RCi*SI  N  (  AM(  ) 

110  CALL  PLOT (XX, YY.2) 

DO  112  1=1,6 
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ANG= ( I-i )*P 1/6. 

ANGS-At'O  ■- 1 
ANGF=  A!JG 

I  F(  I .F0.2*( I  /?  )  )  on  I  )  1 

A  MG 5= A MG 
A  NGF=ANG+3I 

111  COM!  IMUF 
XX=PP*COS(AMC?) 
YY=LP*5IN(AM05) 

CALK  PLOT (XX , Y”  ,3 ) 
XX=L‘'*COS(  AMGF) 

Y  Y=  L  :>*S  I N'  (  AM GF  ) 

112  CALL  PLOT ( XX , Yv ,2 ) 

***  patt-ln  plot  *** 

:.)o  i2o  ip=o,:-<sn,iDHs 

I=IC+ | 

ST'-<-CABS(  HT(  I  )  )/FK'X 
IFCIPLT-2)  121.122,123 

121  2n=L'P*HTM 
00  '10  125 

122  P0=h  P*HTV*rTV 
GO  TO  125 

123  I  F  C  E'l  V .  IT  .0.0 1  )  LTV =0 . 0  I 
kD=?C.*AL00!  0(HTV  ) 

I F  ( t'O  •  LT .  -40.  )  RP=-40. 
20=f;P*(  i<n+40.  )/40. 

125  COM2  I MUF 

A M0=  I  p/I'PP 

xx-h:  -*sin  c  Ar:c ) 

YY=LO*COS(AtlO 
I PE*- -p 

I  F( I .HO.  I  )  I PFM=3 
1 2  C  CAL!  PLOT (  XX,  YY,  I  PHM) 

CALL  PLOT (4. 25, -5.5,- 2) 
130  CONTI HUF 

CALL  PLOT (0. ,0. ,999) 
RFTPHN 


